Ab initio charge, spin and orbital energy scales in LaMnO.3 
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The first-principles SIC-LSD theory is utilized to study electronic, magnetic and orbital phenom- 
ena in LaMnC>3. The correct ground state is found, which is antiferro orbitally ordered with the spin 
magnetic moments antiferromagnetically aligned. Jahn- Teller energies are found to be the largest 
energy scale. In addition it is the Jahn- Teller interaction which is the dominant effect in realising 
orbital order, and the electronic effects alone do not suffice. 

PACS numbers: PACS 71.27. +a, 71.28.+d 



There are several transition metal compounds in which 
the orbital degeneracy is broken spontaneously. Exam- 
ples are KCuF3,0 V203,@ and the manganites,0 which 
are the subject of the present study. In the manganites, 
the crystal field associated with Mn06 octahedra splits 
the manganese d levels into a lower lying t2 S triplet and 
an upper e g doublet. The t^g states are highly localized 
whereas an electron in one of the e 9 states is potentially 
itinerant. The Mn 4+ ion in CaMn03 has a fully occu- 
pied majority t2 S manifold and empty e g states, which 
form a strongly localized core spin S=|. The Mn 3+ ion 
in LaMn03, on the other hand, has an additional d elec- 
tron which, due to the strong intra-atomic exchange, pop- 
ulates one of the e g states, forming an S — 2 spin, and 
which gives rise to a Jahn- Teller (JT) instability. In this 
system each of the oxygen ions is a neighbour to two Mn 
ions and hence the local distortions of the lattice must 
be arranged in such a way as to minimize the energy for 
the whole crystal. This gives rise to Mn-0 bond lengths 
of 1.90 and 2.18A within the manganese oxygen plane, 
compared with 1.96A for the hypothetical cubic system. 




The e g states rotate to form an orbitally ordered lattice 
of d^ x 2_^ r 2 and d^ y 2_ r 2 orbitals in the manganese oxygen 
plane0, shown schematically in Fig. ^ 

This paper reports a first-principles study of charge, 
spin and orbital energy scales in LaMn03, with the 
emphasis on orbital ordering (OO), based upon the 
self-interaction corrected (SIC) local spin-density (LSD) 
theory, which allows d-electron localization to be dis- 
tinguished from itinerancy, in an ab initio manner. Q 
Not only does the SIC allow the localization of an orbital 
of any symmetry, but the SIC-LSD total energies can 
be minimized with respect to the number of localized 
orbitals and their orientation. In addition, minimization 
with respect to the number of localized orbitals yields va- 
lency which is defined as an integer number of electrons 
available for band formation, N va i = Z — N core — Nsw> 
where Z is the atomic number, N cor . e is the number of 
core (and semi-core) states, and Nsic the number of self- 
interaction corrected (localized) states. Furthermore, in 
the SIC-LSD the occupation of each localized orbital is 



FIG. 1: Simplified view of the manganese oxygen plane in 
LaMnOa. The small black circles represent the manganese 
atoms, while the open circles mark the positions of the oxy- 
gens. The distortion has been exaggerated to show clearly 
that the manganese oxygen bonds alternate between long and 
short in the x and y directions (designated Mn(x) and Mn(y) 
respectively), which is associated with the antiferro OO of the 
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determined self-consistently. When localizing an e g or- 
bital its symmetry is automatically broken leading to a 
nonzero value for the local orbital order. For example if 
the orbital that is localized corresponds to d^ z 2_ r 2 then 
the occupation of this orbital and its partner, d x 2_ y 2 , will 
be different. This difference, because of hybridization ef- 
fects, will not be equal to unity as it would be in a fully 
localized picture. 

The SIC-LSD provides the quantitative comparison of 
all energies for LaMn03. This material has a distorted 
cubic structure with OO as shown in Fig. ^and is an A- 
type antiferromagnet (A-AFM) such that the moments 
in the x-y planes shown in Fig. ^are ordered ferromag- 
netically and the planes are stacked antiferromagnetically 
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TABLE I: Summary of all energy scales in LaMnC>3. 



Quantity 


Energy gain 
(mRy per formula unit) 


Electronic energy gain 
due to the distortion 


~ 40 


Localisation energy 


~ 20 


Purely electronic contribution to 
the orbital ordering energy 


~ 5 


Relative stability of 

the A type antiferromagnetic 

state compared with ferromagnetic 


~ 1 



along the c axis. Although the dominant distortion about 
a Mn site is tetragonal there are other components so 
that the net symmetry is lower. We have established 
that the relevant energies are: change in the total energy 
due to the distortion, localization energy, the purely elec- 
tronic contribution to the orbital ordering energy (i.e., 
the change in the electronic energy induced by 00 in 
the absence of distortion ), and the magnetic ordering 
energy; these energy scales are summarized in Tabic H] 

The largest electronic energy is that due to the dis- 
tortion i.e., the gain in energy by imposing 00 in the 
distorted crystal structure. The second largest energy in 
the problem is the localization energy which is the dif- 
ference between the energies of localized and delocalized 
e g states. The value is about 20 mRy per formula unit, 
depending slightly on the state of magnetisation and dis- 
tortion. The magnetisation comes last, as the smallest 
energy scale, behind the purely electronic contribution 
to the orbital ordering energy. The oxygen displacements 
around site Mn(x) are such that the site symmetry has a 
large, 2e xx — e yy — e zz , distortion, as well as other smaller 
distortions that break the symmetry between y and z. 
The Mn(y) site has a similar distortion. Table [D] con- 
tains all SIC-LSD energies when a given e g orbital has 
been localized on a given site, for A-AFM spin ordering. 
All calculations were done for the experimental distorted 
structure. 

The following points may be seen from this table. 
First, the most favourable state is the antiferro orbital 
ordering (AF 00) on the correct site coupled with anti- 
ferromagnetic spin ordering as observed experimentally. 
Second, the results for scenarios 6 and 7 should be identi- 
cal from symmetry. They were calculated independently 
and the difference between them represents an estimate 
of the error on our calculations; the same is true for sce- 
narios 9 and 10. Third, the size of the total electronic 
energy favouring the orbital ordering is found from the 
most and least favorable AF scenarios 1 and 4, which dif- 
fer by ~ 78 mRy, half of which is then the energy gain 
due to 00, i.e. the largest energy scale of Table |U More 
explicitly, the states in scenarios 1 and 3 have identical 



TABLE II: Results for several scenarios where the three t2 9 
orbitals and an additonal e 9 state are localized on the man- 
ganese atoms in fully distorted A-AFM structure of LaMn03 . 
The first column gives the numbering of the scenarios which 
involve AF OO and ferro (F) OO. Columns 2 and 3 indicate 
which e 9 orbital is localized on the manganese atoms with 
the long manganese- oxygen bonds in the x and y directions 
respectively. The relative energies AE, with respect to the 
ground state with the localized d 3x 2_ r 2 orbital on the Mn(x) 
sites and the d 3y 2_ r 2 state on the Mn(y) atoms, is displayed in 
column 4 (A-AFM). Column 5 refers to relative energies AEd 
obtained from a localized model which incorporates only the 
orbital ordering in a purely tetragonal distortion, (see text) 











A-AFM 




# 


Mn(x) 


Mn(y) 


AE (mRy) 


AE d {mRy) 


o 


1 


d 3x 2_ r 2 


d 3y 2_ r 2 








o 


2 


d z 2 —x 2 


dy2 — Z 2 


8.5 


20 


< 


3 


d 3y 2_ r 2 


d 3x 2_ r 2 


65.2 


59 




4 


dy 2 — z 2 


d z 2 —x 2 


78.1 


78 




5 


d x 2_y2 


^x 2 — y 2 


27.2 


20 


o 


6 


d 3x 2_ r 2 


d 3x 2_ r 2 


39.7 


29 


o 


7 


d 3y 2_ r 2 


d 3y 2_ r 2 


40.5 


29 




8 


d 3z 2_ r 2 


d Zz 2_ r 2 


49.5 


59 




9 




^z 2 —X 2 


56.6 


49 




10 


dy2 — z 2 


dy2 — Z 2 


57.3 


49 



pattern of antiferro orbits. The difference between them 
is that in the former case the lobes of the orbits match 
the lattice distortion, while in the latter case the orbit 
is misaligned with respect to the distortion. Band struc- 
ture effects favour these states equally and the difference 
in their energies represents the importance of the lattice 
distortion for the 00 energy. 

The energy scale that characterises the stability of A- 
AFM versus FM is the smallest energy scale, at ~ 1 mRy, 
which is similar to energies that gave correct exchange 
constants for NiO(lOO) surface. |lOj Thus, the ordering 
sequence of the rows will be the same for the A-AFM and 
the FM structures, which indicates that the dominant 
cause of the orbital ordering is not the A-AFM magnetic 
structure, as has been postulated. [ill ll2T] 

It is instructive to see how far a localized model that 
assumes only a locally pure tetragonal distortion can ac- 
count for the ab initio results. The energy per site is 
given by u — uq — U£>cos29 where uq and ud are the en- 
ergies that do not and do depend on the distortion respec- 
tively. The angle 9 defines the orbit. For the site Mn(x) 
we have 9 — : d^ x 2_ r 2, 9 — ±7r/6 : d x 2_ y 2,d x 2_ z 2, 
9 = ±27r/3 : d^y 2 _ r 2 , d% z 2_ r 2 , 9 — tt/2 : dy2_ z 2, and 
equivalently for the Mn(y) site. The energies of scenar- 
ios 1 and 4 may be used to fix uq = up = 39 mRy for 
the A-AFM phase. The results given in column 5 of Ta- 
bleware found using this expression. It is seen that the 
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order of energies obtained by this simple model repro- 
duces the trends seen in the ab initio results. Moreover, 
most ab initio energies lie within 10 mRy of this very 
simple model. It appears that the discrepancy between 
the model and the first principles results is mostly due 
to other distortions than tetragonal. In order to investi- 
gate this we have performed ab initio calculations for a 
crystal with a pure tetragonal distortion (of the observed 
magnitude). In this case we found a slightly increased 
value for ud ( 42 mRy) and the deviations between the 
model and the first principles calculation are reduced to 5 
mRy, for all scenarios. The size of these remnant energy 
fluctuations designates the scale of the band electrons' 
contribution to the total energy. The imposed 00, lat- 
tice distortion and magnetic structure induce changes to 
the conduction electron states, which leads to changes 
of the order of 5 mRy to the total energy. This rep- 
resents the third energy scale of Table HI which is the 
second smallest of the energy scales considered. Another 
estimate for the band-electronic effects was obtained by 
evaluating the 00 total energy variations in the cubic 
phase. Our calculations gave 3.7 and 8.3 mRy in this 
case, for FM and A-AFM ordering respectively, i.e. 00 
is drastically suppressed without the lattice distortion. 

To proceed with the analysis of the energetics of 00 
in the manganites, we now consider the lattice distortion 
mode. The strength of the orbital-lattice interaction may 
be deduced from the measured values of the elastic con- 
stants and the size of the lattice distortion. To see this, 
we follow Ref. and consider an expansion of the to- 
tal energy of an electronic system at T =0 in terms of 
a local distortion in the unit cell at Ri, e^.The restoring 
force for that mode is characterized by the force constant 
K, and the elastic energy of the distortion has a simple 
harmonic dependence on amplitude. Countering this is 
the nonzero orbital order, which contributes a negative 
energy: 

Utot = U + U e i{{e,}) + ^Y, K ^f 

i 

i i 

In deriving this expression we assumed that only linear 
terms in ei were required for the 00 energy, and defined 

dUel({ti}) 

de„ 

At equilibrium, the elastic and 00 energies must balance, 
and the distortion parameter is found by minimising the 
energy, i.e., e_ = Therefore the total lowering of 
energy due to the lattice distortion and 00 is found by 
substituting the above expression for e_ into equation IjTjl: 

AU = U e i(0) -\Y, = U * (°) - 5 _C Ui6i - 

i i 
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FIG. 2: The orbital and spin resolved densities of states cor- 
responding to scenario 1 (see Table IHl . i.e. with the e 3 lobes 
oriented along the long Mn-0 bonds. The top panel shows 
the majority spin partial p densities of states for the O(x) 
sites, while the middle and lower panels show the majority 
and minority (summed over all d orbitals) spin partial d den- 
sities of states for the Mn(x) site respectively. The energies 
(in Ry) are relative to the Fermi energy at E = 0. The partial 
densities of states are in units of states/spin/Ry. 



The distortions are known in the manganites, e_ ~ 0.30 
A, j_| as are also the optic phonon frequencies, leading to 
K ~ 12.5 eV/A 2 .[]]| Hence, we can determine an ex- 
perimental value of Ui — 3.5 eV/A. The SIC-LSD theory 
leads to the estimate Ui ~ un/ti ~ 1.8 eV/A, which is 
to be considered in excellent agreement, given the gen- 
eral uncertainty, not least in estimating the appropriate 
value of -K".[___J Otherwise stated, neglecting U e i(0), a to- 
tal Jahn- Teller ordering energy of ~ \ u d ^20 mRy is 
arrived at, which is in rather good agreement with the 
estimate of 0.48 to 0.58 eV =35 to 43 mRy of Ref. [U 

To shed more light on the hybridization effects involved 
in 00, we inspect in detail the magnetic m and orbital 
I quantum numbers resolved densities of states for two 
00 configurations, namely scenarios 1 and 3 of Table ITT1 
which are shown in Figs. El and [3 respectively. Localised 
orbitals are seen to have only marginal admixture (below 
0.1 electron) of the delocalized e g state. The latter (in 
green) can be found both in the valence and conduction 
bands. The localized orbitals do still contribute some 
weight in the valence and conduction bands of the Mn(x) 
majority DOS in Figs. [2 and [3 In particular, t 2g states 
contribute to the valence and conduction bands. This 
reflects the slight hybridization of the t2 g states, show- 
ing that they are not of pure ionic character but acquire 



4 





- Py 

- P z 

- Px 

— i — f " A~ 


-3/-r 2 

11 




— minority spin d 





-0.8 -0.4 0.4 



FIG. 3: The same as Fig. Hbut for scenario 3 (see Table |TTJ, 
i.e., the e 9 lobes are oriented along the short Mn-0 bonds. 
The top panel shows the majority spin partial p densities of 
states for the O(y) sites. 

some covalency. For the localized e g state this hybridiza- 
tion increases and is most pronounced for d 3y 2_ r 2 in Fig. 
|S| Here we find d 3y 2_ r 2 states in both the conduction 
band and at the top of the valence band, summing up 
to 0.2 electrons. Obviously, this localized state, whose 
lobes point along the short Mn-0 bond length, hybridizes 
strongly with the p x orbital of O(y) site and this can be 
seen as the p x weight in the conduction band of O(y) 
partial DOS in Fig. |31 In contrast the d 3x 2_ r 2 weight in 
the valence and conduction bands for scenario 1 is less 
than the d 3y 2_ r 2 weight in scenario 3. In scenario 1 the 
e ff orbital is more localized: 0.82 localized d 3x 2_ r 2 elec- 
trons versus 0.74 d 3y 2_ r 2 in the DOS peak around -0.8 

By. 

These results can be easily understood in terms of the 
OO reducing the overlap between the localized (SIC) e g 
states and the 2p states of the adjacent oxygen atoms. 
This is reflected in the DOS of scenario 3 by the increased 
hybridization in the O 2p channel around -0.8 Ry, over 
the whole width of the valence band and even in the con- 
duction band (Fig. in comparison with the DOS of 
scenario 1. Thus localized e g states which have lobes 
along the long manganese oxygen bonds will be energet- 
ically favourable configurations, while those that point 
along the short manganese oxygen bonds will correspond 
to energetically unfavourable configurations. Therefore 
the next two in the sequence of favourable scenarios are 
also antiferro OO involving in one case d z 2_ x 2 and d y 2_ z 2 



(8.5 mRy above the ground state) and in the other case 
the d 3x 2_ r 2 and d 3y 2_ r 2 but with ferromagnetic spin ar- 
rangement (unfavourable by ~ 1 mRy). It is significant 
that even with the ferromagnetic spin arrangement the 
antiferro OO is lowest in energy indicating again the im- 
portance of the JT effect. 

In summary, using the SIC-LSD theory it has been pos- 
sible to investigate the orbital, spin and charge ordering 
of distorted LaMn03. We find that our calculated val- 
ues for the orbital energy depend strongly on the lattice 
distortion and are essentially independent of the mag- 
netic order. We have used various numerical estimates 
to support our claim that the Jahn Teller interaction is 
the dominant effect in producing OO. The lattice effects 
are big enough to account for the observed OO, meaning 
that it is not necessary to invoke additional contribu- 
tions of an electronic origin. This, however, should be 
compared with the results of other papers, [lil|l2j where 
it is claimed that OO could occur from electronic effects 
alone. We agree with this in so far that we also find 
an electronic effect of the correct symmetry. However, 
we find that the energy associated with the pure elec- 
tronic effect is relatively small and that the size of the ef- 
fect associated with the distortion is large enough to give 
rise to the structural transition. In addition, this could 
be corroborated by comparing with independent experi- 
ments. This is an important finding for the transport in 
the doped materials because as the electrons become de- 
localised the lattice is unable to respond fast enough and 
one reaches the large polaron regime. On the other hand 
OO that is purely electronic could coexist in the metallic 
phase giving rise to residual OO and extra contributions 
to the scattering. 
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